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Die Anwendung der elliptischen Functionen bei Berechnung 

absoluter Störungen. 



§1- 

m 

Die hier gezeigte Verwendung der elliptischen Functionen in der Störungstheorie 
ist nach dem Grundgedanken bereits im Jahre 1886 von mii' zum ersten Male auf- 
gestellt und damals auch schon Herrn Weiei'Strass mitgetheilt worden. Im Jahre 1889 
habe ich dann die sechs Grundformeln im Bulletin Astronomique bekannt gegeben. 
Ich wende bei Einfuhrung der Formeln die mir geläufigeren Weierstrass'schen Bezeich- 
nungen an. Demnach sind es folgende Relationen, von denen ich ausgehe: 

/(Js 
Vi (a—ci) (s—e-j) (s—e-O 

pu 

(ei — 6-2)0 -u + Ci^u = O'i^u 

(et — ei)a^u -\- Oi^u = afy 

(ea — ei)ahi + a^hi = Oiht 

(e'2—ei)ai^u + (ci— e^aa^w = (ei—ei)a2^u. 

Wenn nun femer unter e die excentrische Anomalie eines Planeten verstanden wird, 
so setze ich 

sm-^« = \ei — 62 '-—y cos-^« = -*-, 1. 

wo die zusammengehörigen Werthe e und u sich durch die Bedingung bestimmen, 
dass gleichzeitig 

e = 0, u = 0. 

■ 

Aus den Gleichungen 1) folgt nun sofort 



au • <r.u 



Sin« = 2Ve,-c» • -^^ 
= 24?-- 1 



cos« 

a^'H 

au 



1 a^u de ^ — — — d 



2. 



— 2 — 

Nun beweist man leicht die Richtigkeit folgender Fonneln: 

d au a^ua^u 

du a^u a^ua^u 

d a^u . V aua^u 

niit deren Hilfe man die erste und dritte der Gleichungen 2) so reducirt: 

• 

2V«i — «a d a^u 



sm fi = 



e^ — «3 du a^u 



du — ^^^^ ^* a.u 



Der Radiusvector des Planeten ist gegeben durch 



— = 1 — ecos€. 

a 



Transformiren wir durch die Variabele ti, so kommt 



a \ 1 + « «^a^/ 



Auf Grund dieser Gleichung nun bestimmen wir das eingeführte System elliptischer 
Functionen näher so, dass es- mit dem astronomischen Problem in engeren Zusammen- 
hang tritt. Es besteht nämlich auch die Relation 



«1— «8 "a''^ ^1.-^8 V^ 

und die Grössen ei, es, es genügen bekanntlich iöimer der Gleichung 

ei + 62 + ez = 0. 

Setze ich daher jetzt 

ei— 63 = 1 

62 — ^3 2e 9- 

€1 — 63 ~ 1 + 6' 

so sind also im Ganzen drei Gleichungen vorhanden, um die drei Grössen ei, «2, «3 
als Functionen der Excentricität des Planeten zu bestimmen. Zur Abkürzung wird 
noi^h gesetzt: 

62—63 ^ 2 

61 — e-i 



— 3 — 

m 

Dann bestehen folgende später häufig zu brauchende Foi'meln: 

2e „ 1 — e ^ 
— * > i = * 



1 + e ' 1 + 6 



X2 



e = 



1 + e = -r~:r, 3^ 



1 + x'2 ' ^ 1 + x'« 

2x' . 2x'3 



Vl-e« = :i^^-,T, 1-e = 



1 + X'» ^ 1 + x'*' ' 

Setzen wir femer noch zur Abkürzung 



ffgW 



= XU, 



so stellen sich unsere vorhin gegebenen Formeln leicht, wie folgt, dar: 

dB ^ 2Wxu\du 

2x' dx 1 + x« 2»'* « 

* 

Endlich bemerken wir noch, dass, wenn g die mittlere Anomalie bedeutet, 

dg = ^dB, 
wonach also 

^9 = f:p?-»x'(«)-^w- **• 



§ 2. 

Die Integration der Störungsgleichungen in der Form, wie ich sie. nachher 
anwende, erfordert noch die vorgängige Berechnung der Ausdrücke 

rcos/, rsinf; -^ , -^; sm/, cos/+e, 

wo / die wahre Anomalie des gestörten Planeten, r dessen Radiusvector und e seine 
Excentricität bedeutet. Wenn man sich nun der Formeln erinnert 

~ sin / s= Vi — e^ • sin € , ~ cos f = cos e — e , 

SO erhält man sofort 

a ' 1— X* du 



r - 4«'* 2*'« ,, V 

— COS /^ = -z jr- 5— X^Ou) 



und aus diesen Formeln in Verbindung mit denjenigen des vorigen § findet sich 
weiter : 

« sinf __ 1 + x« 1 rf/ 

acosf _ 1 + x'i JL (l + /y J^ 

^ 6. 

. , 2 1 dx 

1 — x'* x^** ^^ 

r X 2 1 4x'* 

Dies ist aber noch nicht die endgültige Form, in der wir diese Grössen benutzen. 
Denn die Störungsgleichungen erscheinen alle in der Form 

WO F von den Coordinaten des gestörten Planeten und den partiellen Ableitungen 
der Störungsfunction abhängt, und h die aus Hansen's Theorie bekannte Grösse ist 

, an 

h == 



cos (p 



{n = mittlere tägliche Bewegung, tp = arc sin e). Es tritt nun aber u an Stelle von t 
als unabhängige Variabele; und es ist 

d^ d^ dt _^d^ 1 dg 

du dt. du ' dt n du 

Wir müssen also alle unsere Ausdrücke ^- ^ mit ^ multipliciren. Dies geschieht, 

indem wir die als Factoren auftretenden Grössen 6. mit j- multipliciren. Diese werden 

daher jetzt als rationale Functionen von %u und /w erscheinen. Um dieselben später 
berechnen zu können, führe ich einen Satz an, den ich schon 1889 im Bulletin 
astronomique gegeben habe. 



§ 3- 

Unsere Function xu genügt der Differentialgleichung 

Hieraus folgt durch Differentiation, nach Weglassung des Factors -p- 

g = (1 + x'Oz - 2*h\ n. 
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Machen wir in I. die Sabstitution 

dx ^ ±_ d_x* 

du 2x ^** 

und diflFerentiiren, so kommt, wieder nach Weglassiing des Factors -j^ , 

gv = - 2 + 4 (1 + Ox' - 6*'V. m. 

Differentiiren wir II. und III. zweimal und reduciren, so erhalten wir 

= ((1 + X") + 12 x")x - 2üx" Cl +.x'^) X' + 2ix'*x' IV. 

^ = — 8 (1 + »'■)+ (16 (1 + X'«) + 72x'») z' — ISüx'^Cl + x'»)x* + 12üx'*x*- V. 

Die Betrachtung der Gleichungen II. bis V. legt die Vermuthung nahe, dass 
allgemein i.st 

~f , VI. 



,, = ^^->X + Äi(-^X' + . . . + ^n<''^X""+ * 



rfw 



ün 



= Su("> + Ä<"^X* + . . . + Bn^^^X^*"^^' 



In der That lässt sich die Richtigkeit dieser Gleichungen durch den Schluss von n 
auf w + 1 nachweisen, wobei sich dann noch die Recursionsformeln ergeben: 

j(- + i) ^ (14.x'2)(2r + l)Ut"^— x'2.2r(2r— l)41i — (2r+2)(2r + 3)^i"ji 
S(n+i) ^ 4(l+x'2)r»5<-^— x'2.(2r— 2)(2r— l)ßj.li — (2r + l)(2r + 2)B^-|i 

Die durch VI. gegebenen Sätze lassen sich leicht in Worten ausdrücken; mit ihrer 
Hülfe ist es also möglich, die geraden und ungeraden Potenzen von xu aus xu und 
X^n und deren Ableitungen zu berechnen. Ganz dasselbe gilt für diene gativen 
ganzen Potenzen von x^* vennöge der Relation 

—1 

3C(m) = x'X (m + w) , 

wo 

00 






ß—ez) 
«i 
Nun ist 

-.x'x« = 1 + ^^.(-irr^cos— , 

wo ^^^ die bekannte in der Theorie der elliptischen Functionen auftretende Grösse ist. 
Die Ent Wickelung für x*w erhalten wir aus der Relation 

d^ 



— 6 — 
und zwar wie folgt; 

\n / 1 g^r l g9r ^ 

Den Modul q berechnet man bekanntlich aus 

WO 

1 + \*' 

Setzen wir noch 



' *=7r"''' ^ = (jy(i + x-) = |(«' + .?'). « = 



so können wir folgende Formeln für die ersten sechs Potenzen von %u aufstellen: 

{&iu) =00+ 4 ^ ( — l)*" Qr COS rv 
(«X«)* = 6o — 85 ( — l)' hr COS rv 
(«X«)* = Co + 16 •^ ( — l)*" Cr COS rv 
(txu)* = do— S22 (—lyUr COS rv 
(ejtu)» =eo+ 64 5 ( — l/ er cos r» 
(«X«)« = /o — 128 5 (— 1)' fr cos r» , 

00 = 1 ; 5o = l-85(-iyj-^; co = 2A 

do«|-^6o-i«?*; eo = I »?« (l-x» O ; fo = j^do~eW(h-^^} 

dr==(^l + f')br ; «r={^eo + r«(AH-r«)}«,;/, = irf.-l(,J-f5r«)ft. 

Die Grössen €, i; sind gegeben durch die Reihen 

e = (1 — 2q + 2q^ — 2q^ + — . . .y 
f] = (l+2q + 2q^ + 2q'+ )\ 

die bei den beabsichtigten numerischen Anwendungen in unserem Planetensystem 
stets mit dem dritten Glied abgebrochen werden können. Mit Hülfe der obigen 






• 'V 



Darstellnngen der Potenzen der Function x^ ist übrigens gleichzeitig auch eine solche 
für die Potenzen von xu erlangt, bezw. für alle Ableitungen von ^a. In der That 
kann x w selber aus der Reihe für x^ stets durch directe Differentiation abgeleitet 
werden; und für eine beliebige ganze Zahl n ist 



(: 



^) " = (-1)- (1-1 + x'^ • X' + ^''x'y , 



d. h. diese Potenz von x'^ erscheint als ganze rationale Function von x^^ 

wo die Qr Functionen des Moduls x' sind. Und aus dieser Formel fliesst dann sofort 
auch für die ungeraden Potenzen von x'u folgendes Ergebniss ; 

(X'«)««+» = (-1)- 1 (z + y a,z^ + y «»Z-^ + . . . + -^^ a,„x*-+i) , 

Des Zusammenhangs wegen sei noch bemerkt, obgleich im Folgenden kein 
Gebrauch davon gemacht werden wird, dass man in gleich einfacher Weise auch für 
die negativen Potenzen der Function xw ganz analoge Sätze ableiten kann, wie sie 
vorhin für die ganzen positiven Potenzen erhalten worden sind. Man wird dazu, ' 
indem man setzt 



X(n) = ^, 



ausgehen von den Gleichungen 






rfw» 



= 2x'>+4(l +0^ — 6|\ 



und mit diesen dann ebenso weiter verfahren , wie es oben mit. den entsprechenden 
Gleichungen für %a geschehen ist. 



Nach diesen Auseinandersetzungen über die Darstellung der Potenzen von 
XU und x'm ist es nun möglich, näher in die Anwendung einzutreten. Zunächst 
können nun die Formeln 5. und 6. in § 2 in der dort erwähnten Weise vervollkommnet 
werden. 
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Wir finden 



r , j, dg 4x' <^ • / n* 

T^^'^'^-Sr = »'(i + x-y w^"'^«) 

a . M dg 4x' 1 , 

r du X* /m X* '^ 

Die dritte dieser Gleichungen scheint die Einführung der Function log^w zu ver- 
langen, da sie die logarithmische Derivirte von xu enthält. In Wirklichkeit .ist die 

Anwendung dieser Function aber nicht erforderlich. Denn — sin / ist in den 

T 

Störungsgleichungen multiplicirt in r • -^^ , wo ß die StÖningsfunction bedeutet; und 
diese letztere Grösse werden wir darstellen als ein Aggregat homogener Functionen 

TT 

von — sin/ und — cos/, vom zweiten und den höheren Graden, so dass also bei 
a a 

Herstellung der definitiven Form der Störungsformeln die negative erste Potenz von 

XU herausfällt. 

Bevor diese letzteren Formeln gegeben werden, möge noch darauf hingewiesen 

sein, dass ich in den Rechnungen über absolute Planetenstörungen besondere Ent- 

Wicklungen für die Componenten der störenden Kräfte, also -=w > ** -^ > wenn man sich 

zunächst auf die Störungen in der Bahn beschränkt, nicht anwende, sondern diese 
beiden Ausdrücke aus der Entwicklung der Störungsfunction wie folgt berechne. 
Diese letztere ist natürlich herzustellen als eine trigonometrische Reihe nach dem 
Argumente w, sowie anderen Argumenten, welche irgendwie von dem störenden 
Planeten abhängen; die Integration der Störungsgleichungen ist nach der Variabelen « 
auszuführen, welche an Stelle der Zeit t jetzt als die unabhängige Veränderliche des 
ganzen Problems erscheint. Man wird demnach haben: 

?fl ?fl df . dii dT 



= ^^ 4. 



du df du Dr du 

Die' hier nothwendigen Ableitungen von / und r ergeben sich leicht, die 
letztere ohne weiteres aus Gleichung 4. §1, die erstere durch Beachtung der Relation 



df= :i^.n-e''dff. 
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Es wii'd erhalten 

-^ = A *!1 = 2r.— 
du %u ^ du xu^ 

wonach sich dann folgende Gleichung für die Transformatiop von -^ in -^ ergiebt 

Hier wird man also nnn -x— durch directe Düferentiation aus dem entwickelten ß 

cu 

erhalten; während andererseits sich auch zeigen lässt, wie mit Hülfe einer ganz 

einfachen Anwendung des Euler'schen Satzes über homogene Functionen auch r -^ 

aus der Entwicklung von ß ohne jede Mühe hergeleitet werden kann. Dieser Nach- 
weis wird im § 6 gegeben werden. 



§ 6. 

Die Störungsgleichungen*) lauten nun, wenn wir uns hier zunächst auf die 
Berechnung der Störungen erster Ordnung in der Bahn beschränken, also von den 
Breitenstörungen absehen : 

dw „,Öß 

dl orSin/ aß , 2ä ( .f^ , .cos/ , ,„. ^ . A ?ß 

di r dr x^(l + m)\ r / ?/ 

d^ 



dt 

d^ 
dt 



=i w + S* r cos/ + rj' r sin/ 
= ih (S sin/— r](cosf + eU . 



Hier bedeutet x^ die Gauss'sche Constante, m die Masse des gestörten Körpers und h 
ist wieder definirt durch 

h = , 

cos^ 

wo a die mittlere Entfernung des gestörten Körpers von der Sonne, n seine mittlere 
tägliche Bewegung und g> = arc sin e ist. 



*) Ich leite die hier gegebene Form dieser Gleichungen an anderer Stelle ab. 

2 
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Nach Auflösung der Störungsgleichungen berechnen sich die wahren Coordi- 
naten r, v des gestörten Köi'pers in seiner Bahn aus folgendem Gleichungssystem, 
in dem noch Mq die mittlere Anomalie und ^o die Perihellänge der Epoche bedeutet: 

€ — c sin € = wC 4- ^Wq 
r C0S2; = a cos€ — ae 
r sin t; = a cos qp sin e 

= i? + ^0 
X = r(l + v). 

Die drei ersten Störungsgleichungen sind auf unsere unabhängige Variabele a 
zu transformiren, wonach sie folgende Gestalt erhalten, wenn wir noch zur Abkürzung 
der Schreibweise 



setzeö. 



dw 
du 



x^»x^2 dl 

1 + X'2 du 



er du 






Behufs Integration dieser Gleichungen sind nun P und Q in trigon. Reihen 
nach den Yielfachen des Argumentes u zu entwickeln. Zu dem Zwecke betrachten wir 
zunächst die Störungsfunction ß selber, die in bekannten Bezeichnungen die Form hat 

ß = I? i-j — -^ cos (riryj . 

Der wichtigste und für die la,ufende astronomische Arbeit am meisten in 
Betracht kommende Fall im Sonnensystem ist derjenige der Störungen der kleinen 
Planeten bezw. der periodischen Cometen durch Jupiter und Saturn, mit Rücksicht 
auf welchen Fall auch vorliegende Methode in erster Linie zum Vorschlag gebracht 
wird. Ich setze deshalb auch voraus, dass r>r. Dann lässt sich für ß folgende 
bekannte Darstellung geben: 
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^ m' (1 r^ r^ ^ \ 

Indessen ist zu beachten, dass wir mir die partiellen Ableitungen von fl nach 
den Coordinaten des gestörten Körpers brauchen werden und dass wir uns in Folge 
dessen bei den aufzustellenden Entwicklungen auf die Betrachtung des Ausdruckes 

iii = H-~i~z: i~^^2 H — 77X3 + • • •/ 

1 + «I \r^ r* / 

beschränken dürfen. 

Die Kugelfunctionen Xi , X-.. , ... sind definirt durch die Gleichungen 

Xi = Xi = cos (r^r') 

A2 — -g-Ai — y 

X3 = "o""^! — ^2 ^ 

A4 — 27^ A, — 2 -^^Ai -I-2T4 

und die Recursionsformel 

(n + l)X«+i-(2a + l)XiX« +nXn-i = 0. 

Wir eiünneni nun, dass mit i, i' die Neigungen des gestörten bezw. störenden 
Körpers gegen die Ekliptik (oder die sonst gewählte Fundamentalebene), mit *, ^' 
die Längen der aufsteigenden Knoten der beiden Bahnen in Bezug auf dieselbe Grund- 
ebene bezeichnet werden, und dass 

CO = 7t—» ,' cy' = n—»'. 

Ferner sei J die gegenseitige Neigung zwischen der Balin des gestörten und der- 
jenigen des störenden Körpers. Dann definiren wir noch zwei Grössen ?>, t//, deren 
geometrische Bedeutung aus den Definitionsgleichungen erhellt, nämlich 

sin|-J^sin|(()P + V;) = sinlc^-^O sin|(i + 
sin -^ J^ cos y (y -h V') = cos-i-(d— *') sini(i — 
cos 2" J sin — (g> — rp) = sin ^ G'> — ^0 cos yCf + i') 

cos-^J^cosy (9 — 1//) = cos -(*—*') cos^(e— 0; 

und setzen endlich 

2C^= (o — w'—((p—xp) 
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Dann können wir cos (r, r') = Xi in der Form schreiben : 

Xi = cos3ljcos(/-/ + 2f/) + sin2|j^cos(/+/+2r) 

und hierfür lässt sich endlich, wenn man 

^\ = cos^i-J^cosC/ — 26^) + sin^i- Jcos(/ + 27) 

IV = cos^i-Jsin (/ — 2C7) — sin2|-J^sin(/4-2F; 

einführt, nehmen 

Xi = l'i cos/ + l'a sin/, 

wo man noch bemerken wird, dass S\ > I' 2 die Werthe von Xi für /= und/=^7r sind. 
Diesen letzten Ausdruck von Xi führen wir nun in die X„ ein und bezeichnen noch 

— cos/ = X. — sin/ = y . 

Dann haben wir für den hier zu verwendenden Theil der Störungsfunction die Ent- 
wicklung 

a£ii = 82,0 x^ -f Si.i xy + S0.2 »^ « 

4- 0^8,0 a?3 + ^2,1 x'^y Ar Sy% xy^ + Ss,o y"^ 

4- .... =2Sm,na;'^»'* 

oder 

afli = Sl^x, 

wenn //^ eine homogene Function x^^^ Grades (von x = 2 an) der Argumente a?, y 
bedeutet. Es ergiebt sich daher schon hier nach bekannten Sätzen das Bestehen der 
Relation 

ar -^ = ZtcH^ , 
vr 

wo hier nach obigem fl statt fli geschrieben werden durfte. 
Was die S-Functionen angeht, so sehen wir, dass 

^M=:iTp;j-;i\-2"^^ —■ 2;>^*'*==T+^7^'^^'^'' ^«'2 ==1+^7^(2"^^ — 2Jn.s.w., 
wonach also Sm,H allgemein sein wird 

und B,n,H sich bei der Entwicklung von ii ergiebt, wenn in dieser Function r' = a 
gesetzt wird und alle Glieder, in denen m + w <: 2, .weggelassen werden. Diese 



2-« 
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Coefficienten werden daher als Entwicklungscoeflicienten auftreten, wenn wir den 
Ausdruck 

(1 + «» + j^^- 2r.» - 2^2^)"* 

durch eine Potenzreihe nach den Argumenten x, y darstellen (und dabei wieder die 
Glieder unterdrücken, für die m + » <; 2). 

Wie diese Entwicklung anzustellen ist, ist bekannt. Wir erhalten folgenden 
Ausdruck für ©»,,»: 

S (2m + 2»-!)! 1 / (2»t4-2n-3)! 

""'" ~ 2"+»-»(»i + n— l)!»i!w! ^' ^' 2 V 2'»+-»('« + «— 2)! m! (»— 2)! ^' *' 

4. (2>» + 2«-3)! ^ m-ipn\j^L( 2/n + 2>t-5)! ri^r»-* 

1- 2'"+"-« (m + «—2)! (»1—2)! n! ^* ^' / "•" 4 l 2"+—» (tn + »—3)! m! (n— 4)! 2! 

o (2w+2n-5)! ' 1'^»-« ^a"-« (2m + 2n-5)! iiüZ^J^l^ 

"^ 2'»+"-»(m + n— 3)!(m— 2)!(w-2)! 2! "^ 2'»+"-»(m + n— 3)! (m— 4)! «! 2! / 

_ J^ / (2»i + 2«-7)! ri-'r»-" , q (2m + 2» -7)! g'^m-»|>^,-4 

8 \2"'+»-*(m + n— 4)!m!(n— 6)! 3! "^ 2»+— *(m + «— 4)!(m— 2)!(n— 4)! 3! 

. , „ (2m + 2n-7)! ri— * l'2"-» , (2m + 2n-7)! iiül!|'?l\ 

-r » 2"'+"-*(»» + n— 4)!(m— 4)!(n— 2)! 3! ■•" 2'»+"-*(m + «— 4)! (»i— 6)!«! 3!" / 

I • • • • 5 

und wenn wir noch mit a' die halbe grosse Axe der Bahn des störenden Körpers 
bezeichnen, so können wir schreiben: 

Es mögen die Coefficienten für die drei ersten //-Functionen, deren man 
immer benöthigt sein wird, hier kurz zusammengestellt werden: 

©2,0= "2 ^'^^ — "2' ®^'^ ~ 3f'ife', ©0,2 = "2 ^^'^"Y ' 

©3,0 = -2l'l^ "2^'' ®2,1= 2-§'l^r2 -2-|'2, 6l,2 = -2^'l^2^— -g-f'l, ©0,3=-2-^'2^ y?2 

35t, 4 ^^t' 2_L 3 Ä 35 1' St' ^5t/ »., ^ 105 1, 2t' 2 15 t/ 2 15t/ 2_L 3 

4,o = -g 51* — T* "^ 8"' ^3,i = -2 Si §2 — 2-5152, ^2,2 = -;^Si §2' — T^ — 4"52'+-^, 

©' 35 1' t' ^ 15 1/ t' i« 80 t' 4 15 t/ 2 I 3 
l,3=-2'*»2 2'1^25 ^0,4= g^$2 X'^ 's"' 

J J 

§6. 

Die Störungsfunction ist nun also dargestellt in der Form 

aß = 2Sf,,y xf'y^ (fi + y> 2). 
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Die explicite Darstellimg der Prodiicte a;/"«/" macht nach den Darlegungen und Formeln 
der §§ 1 — 4 keinerlei Schwierigkeit. Wir können uns daher zunächst der noch vor- 
liegenden Hauptaufgabe zuwenden, die Coefficienten S als Functionen unserer unab- 
hängigen Veränderlichen u darzustellen. Es ist nun bekannt, dass S^^y sich aus- 
drücken lässt als eine lineare Fimction von cos/*, cos2/', ..., cos (2^ + 2»')/, 
cos (2iu + 2v + 1)/'; sin/', sin 2/', . . . , sin (2fi + 2v)/\ sin (2/£ + 2v + l)/\ Und 
zwar sind die Coefficienten dieser Darstellung ganze rationale Functionen von 

sin ^ -ctJ-, e' und r j;, , wenn e die Excentricität des störenden Körpers bedeutet. 

2 1 — e ^ 

Indessen können uns diese Darstellungen in der Praxis nichts nützen. Als störenden 
Körper haben wir, wie schon erwähnt, vorzugsweise Jupiter und Saturn im Auge. 
Und alle auf diese grossen Planeten bezüglichen Entwicklungen, welche zur Zeit, 
solange umfangreichere Publicationen von Seiten Gylden's und seiner Schüler noch 
ausstehen, allgemein zugänglich sind, sind diejenigen von Leverrier, der als Ent- 
wicklungsvariabele die mittlere Anomalie des störenden Körpers benutzt hat. Wir 
müssen also in Hinsicht auf die Anwendung annehmen, die Functionen S seien als 
Functionen der mittleren Anomalie g des störenden Körpers gegeben, wonach sie 
also als Fourier'sche Reihen in Bezug auf die genannte Variabele zu betrachten sind. 
Die Anführung dieser Entwicklungen an dieser Stelle erscheint nicht geboten; sie 
sind oft gegeben, und über sie etwa Tisserand, Möcanique Celeste I, oder die ver- 
schiedenen classischen Arbeiten Hansen's zu vergleichen. 

Es bietet sicl^ uns also das Problem dar, die Functionen sin xj , cos x/ als 
trigonometi'ische Reihen nach den Vielfachen des Argumentes u zu entwickeln. 

Die früher gefundene Differentialformel für die mittlere Anomalie wird nun, 
wenn wir wieder 



setzen 



u = — , du = - dv 
n rc 






(wo rechts zur Abkürzung noch u geschrieben ist). Und hieraus folgt durch Inte- 
gration unter Beachtung der Bedingung </ = , t; = 

Dabei sei einschaltend bemerkt, dass diese Gleichung, in der Form geschrieben 

^ = 9 + ff/W, 
mit Hülfe der Lagi'ange'schen Reihe dazu führt, den einem gegebenen Werthe g^ also 
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einem gegebenen Zeitpunkte, entsprechenden Werth v zu berechnen. Ferner aber, 
wenn g' die mittlere Anomalie des störenden Körpers bezeichnet, so ist 

dg' = —dg = fidg, 

wo also w, n die bez. mittleren Bewegungen sind. Ist also g'o der Werth von g' für 
den Zeitpunkt, von dem an die Integration beginnt, so hat man 

g' = ^v + g\ 4- ^M («^; (?) 5 

wo if die oben im Ausdrucke für g gegebene Reihe ist. 

Betrachten wir nun hiernach g' als Function von g, so können wir jede 
der hier in Betracht kommenden Functionen f{g) = F{v\ q) = mit Hülfe der 
Mac Laurin'schen Reihe entwickeln 

F{v;q) = F(v;o) + q ^o -\- ^y _l_liZ., ^. . . . 

Was die wirkliche Ausführung einer solchen Entwickelung anbetrifft, so habe ich mich 
überzeugt, dass man sehr schnell und, soweit bei Störungsrechnungeu davon über- 
haupt die Rede sein kann, mühelos zum Ziel kommt, wenn man rein numerisch 
verfährt, also von der Aufstellung analytischer Ausdrücke für die Coefficienten absieht. 
Indessen ist es auch ohne Schwierigkeit möglich, solche analytische Entwicklungen 
durchzuführen, und ich will daher noch die ersten Glieder für sinx^' und cosx^' 
angeben. 

Ich setze x^w = v, dann ist also 

x/ — vv — x/o =r ivg) (y; q) 
und man wird finden 

sin (xg' — vv — x/o) = fi(v;q) 

cos (x/ — vo — xgo) = /2(y;q), 
also 

sin xg' = fi(v; q) cos (vv + xg\) +/2 (y; q) sin (yv + xg'^) 
cosxg' = /a (v;q) COS (w + x/^) +/i (ü; q) sin (yo + xg\) 

Es suid also die Functionen 

/i{v\q) = sin[4>'5p(e;?g)], f^Kp'.q) = cos[4y9)(e; ; g)] 

darzustellen durch die M'Laurin'sche Reihe. Behandeln wir zunächst die erste, so 
sind die Werthe /i(«>(t;; ojg zu ermitteln. Wenn zur Abkürzung gesetzt wird 

w = ^vff {v\q) = ^q>r • sin rv , 
so ist 



t)/i <V/i Dm; 



V«.' 
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wenn wir noch die Ableitungen der Functionen (pr nach q in Lagrange's Manier 
bezeichnen. Für die letzte Formel können wir auch schreiben 

"Y^ = /a {v\ q) • Sq>'r • sin rv 

und dieser gegenüber steht die Formel 

Y^ = — /i {v\ q) • S^pV • sin rv 



Man findet weiter 



3% 

-~y = — /i • 2'fp'r • sin ro + /2 • -2? ^'V • sin rv 



3!/; 



'2 

- = — yi . JS ^"^ . sin rv — f^ • -^qpV • sin rv. 



Schreiben wir die beiden erlangten Ergebnisse so: 

3 sin w 



dq 



Ci cos w 



— :^-2 — = 02 sm w + C2 cos w , 



c^2 
so ist entsprechend 

3 cos t€ ^ . 

— = = — Ci sm u; 

Zq 

32 cos m; 



-^2 = — Ca sin m; + S2 cos «; , 



und wir finden leicht weiter 



3? 



3 



S% Sintis + C3 coSvm; 



3'cosm; •t • • ,0 

— ^-^ — = — ft sm m; 4- 03 cos m; , 
cq^ 

wo 



S 



3S2 ^ 3m? ^ 3Q8 . « 3m; 



■:r- Ci -r— 5 Ca = -7; h S\ 

CO da CO 



'2 3gr ' ^ 3g ^ ' 3g ' 

und es lässt sich leicht durch den Schluss von n auf n + 1 zeigen, dass überhaupt 

3» siuM? 



dq" 



= /Sn sin m; + G, cos w 



3" cos m; ^ . . „ 

— ;r'z — s= —CnSvito + Sn cosm; 

^ _ dSn^i ^ Zw 3CL.1 , ^^ 

3g cq ' 3g ' Sg 
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Beachten wir nun, dass für q = , sin w = , cos m; = 1 , so kommt folgende 
formale Darstellung 

slow = Ci-q + Oi--^ + ...+ Cn-—f + ... 

cos»<; = 1 -I- 5, . </ + 52 • -^- + . . . + 5, • -^ + . . . , (Si = 0) 

wcT die Coefficienten S» , Cx sich als lineare endliche Ausdrücke in sin t; , sin 2t? , ... , 
sin rv darstellen. Es wird sich gleich zeigen, bis zu welchem Werthe von r in einem 
gegebenen S» oder 6» die Grössen sin rv auftreten. 

Zu diesem Zwecke betrachten wir den Coefficienten <pr-i der, abgesehen vom 
Vorzeichen und dem Factor 4v, sich darstellen lässt 

wo 



r 1 + (/^ 



und, da stets <? < 1 , für !/; die Darstellung existirt 



m^i^'f 



»^»^ = ^ 2- [1 + * (?*0] , q}(0) = o. 

Für die Bestimmung von Sn^ On ist die Bildung der Grösse ; nothwendig. Nun ist 

d»»(]p^ = )[d**xp -\- nd**^^tp • C?X + . . . 4" '^(V*X' (^« Kür«e halber lind die Indioes weggrelagflen.) 

Aus der Darstellung von tp erkennt man aber, dass alle Ableitungen dieser 
Function bis zur (r — l)ten Ordnung einschliesslich für g = verschwinden, und es 
sieht sich femer leicht ein, dass man, um den Coefficienten von q^ in den Entwicklungen 
von sin«;, cosm; zu erlangen, diese Functionen nur bis zum Gliede mit dem Index r 
einschliesslich zu differentiiren braucht, so dass also dieser Coefficient sich im All- 
gemeinen als lineare Function von sin i% . . , sin ro darstellen wird. Dabei wird man 
noch beachten, dass für g = ü 



da für g = 



dq^ ^ dq^ ^^^^^ dqr' 

= 1 , X = 1 . 



Für die höheren Difi^erentialquotienten von y\> lässt sich eine independente Darstellung 
geben.*) In der That ist \p als Function von y -= r/ anzusehen. Setzt man 



*) Vgl mein Lehrbach der höheren Analysis, Band I, Differentialrechnung. Berlin 1893. 
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so findet man leicht 

i?x = 2 (-1)' (') «" ('■*-"■) • • ■ ('•«-"-» + 1) «"-'- , 

wo die Summe mit demjenigen Wertlie des Index % abbricht, weicher der in 

T 

enthaltenen grössten ganzen Zahl gleich ist. Unter gehöriger Beachtung dieses 
letzteren Umstandes lässt sich daher für den obigen Differentialquotienten folgende 
formale Darstellung geben : 

WO 



^«=(v)-Kn'-)+®r»")-+- 

und die Glieder wegzulassen sind, in denen negative Potenzen von ^ auftreten. 
Die Grössen — , — x sind bekanntlich definirt durch die Reihen 

^ = (1 + 22 4 2gr* + . . . 4- 2«'' + . . 0» 

"^x = (1 -2g + 2g* + . . . 4- (.-Vf%f' 4- • • •)'• 

Eine independente Darstellung der höheren Difierentialquotienten der aus ihnen 
gebildeten Function X scheint sich nicht geben zu lassen, indessen lässt sich die 
successive Differentiation in ihren Ergebnissen leicht genug übersehen, \vie ich denn 
auch in der wirklichen Ausführung der Rechnung in meinen Arbeiten über die 
Jupiterstörungen einiger kleiner Planeten mit grosser Excentricitat (Typus (33) 
Polyhymnia) bei Bildung der Ableitungen von \p von obiger independenter Dar- 
stellung keinen Gebrauch gemacht habe. 

Die Endformeln, welche man zu verwenden hat, vereinfachen sich auch noch 
beträchtlich, wenn man beachtet, dass i(g) = i( — q)^ dass also die ungeraden 
Ableitungen von X für ^ «= verschwinden werden , und dass V föi' ^ = ebenfalls 
verschwindet. Man hat daher für 2 = (wenn der Kürze halber Differentiale 
geschrieben werden) 

dif = xd^ 

d^ip = xd'^^ 

^^q> = ^d'^xd^f -I- xrf'V 

d^tp = ßrfV^^^ + xd'V^ 

d^g> = bd^xdxp + lOd^xl^xp + xd'tp 

d^q> = Ibd^xd^ + Ibd'^d^tp -f- xd^^ u. s. w. 
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Weiter wie bis zur sechsten Potenz von q wird man wohl nie die Entwicklung zu 
treiben haben. Bei der sehr beti'ächtlichen Excentricität {e = sin 19^ 28' 16") des 
Planeten (33) z. B. wird schon log q'' = 1,8137 — 10. Aus Gründen der Raumerspamiss 
will ich mich hier auf die Angabe der Glieder bis zur dritten Potenz inclusive 
beschränken (was übrigens auch der Führung der Rechnung bei anderen Methoden 
entspricht, bezw. noch etwas weiter geht) und verweise wegen der ferneren Glieder 
auf meine in einiger Zeit erscheinenden Tafeln der Polyhymnia. Man erhält unter 
Berücksichtigung des bisher Gesagten leicht, für // r= , 

dq:' ^ ^ rV^^ dq' 1 + g«'* "^ '''^ "^ dq 1 + qH ' 

wo rechts nach Ausführung der angezeigten Differentiationen noch j « zu setzen ist. 
Wir haben also 

div ^ dtp, 



und ferner 



^q ^^ dq 

'dq^' '^^ dq-' 

^hc ^ d^qi 

Tg-* "^ dq-^ 



smr«; 



sm rv 



•sinr«; 



s, = o, c^ = ^^ 



«2 



Zq ?q^ ' ^ dq^ \ Zq 
Endlich bemerkt man, dass 



^, nZw d'^w ^ Z'%0 /Zw\^ 



(4- r^^) = 1 für r = 1 , und = für jedes r>l 

(d^ ü*' \ 
1-^ ^ r— .1 = 2 für r = 2, und =0 für jedes r>l 
dq^ 1 + ö*V9-o 



i^^ i- — rr) ^0 ^ "~ ^ ^^^' ^'^^ ^^"^ =4-6 für r = 3; aber = für jedes 
q -r q 9 - andere in Betracht kommende r. 

E.S sind also folgende definitiven Ausdrücke zu verwenden: 



*i* 



a 



20 - 



^c; 



-r— = — Izxu • sm t; 
—2 = + 36x/i • sm 2v 



cq 



= 4^ 360x14 • sin 2; — 76 xf/ sin 3t; , 



wo nun wieder %ii für v eingesetzt ist. Wir erhalten so zunächst 

sin m; = — 12xf/ • 5' • sin t; + ISxiw • q^ • sin 2o + (60x^1 4- 216 x^^^) q^ sin v 

_ /^x^u — 72xV)2'sin3t; 

coste; = l-36xV^2 ^ 36x^^2 eos2<; + 108x2u2j3cost;-108xV?'cos3<;, 

welche Resultate wir auch so ordnen wollen 

sin w = Si sin t; + ^2 sin 2,o + «3 sin 3i? ; cos t«; = co + Ci cosü 4-^2COs2ü + Ca cos 3» 

oder 

sin w = ^5i sin Xy, cos m; = 2ci cos i«; . 
Dann ist 

2 sin x/ = ^ (*7 + o) sin [(x^u + X) i; + x^ 0] - -3 (5; — q) sin [(jn^' — X) t; + x/o] 
2 cosx/ = 2(s;i + o) cos[Cxiu + i) t; + x^ü] — 2(5;— Ci) COS[(xiM— X) V + x/o]. 

Damit ist nun aber unsere Hauptaufgabe erledigt. Denn die Coefficienten 5«,„, 
welche oben betrachtet wurden, sind dargestellt durch trigonometrische Reihen nach 
den Argumenten x^/' und können nun also in Functionen der Integrationsvariabein v 
übergeführt werden. Jede einzelne der //-Functionen, aus denen wir die Entwicklung 
der Störungsfunctionen bestehen lassen, also auch ß selbst, stellen wir nun so dar: 

^(^ou; cos [(i 4- i» t; + ce<r] + hf sin [(% + i»t? + /J,.-]} 

und analoge Form erhalten die Differentialquotienten, durch deren Integration die 
Störungen ennittelt werden. Ein näheres Eingehen auf diese Formen, welche nichts 
wesentlich Neues bieten, sondern denen anderer Theorien ganz analog sind und nur 
rechnerisches Interesse haben können, darf hier wohl unterlassen werden. In meinen 
Tafeln der Polyhymnia habe ich alles hierauf Bezügliche in extenso ziLsammengestellt. 
Ich will hier zum Schlüsse nur noch ganz kurz darauf hinweisen, dass die Anwendimg 
dieser Methode mir auch deshalb noch einigermassen vortheilhaft erecheint, weil sie 
mit einer gewissen Leichtigkeit gestattet, solche Glieder, die bei der Integration einen 
sehr kleinen Divisor bekommen, weil fi einem rationalen Bruche nahe gleich ist, 
gesondert zu berechnen. Es ist leicht zu sehen, wie man bei auf diesen Punkt bezüglichen 
Ueberlegungen wieder von den Darstellungen von sin xr/', cos xg' ausgehen und daim ein 
Urtheil gewinnen wird, von welchem Gliede einer solchen Darstellung in der Entwick- 
hmg eine^s der Störungsausdrücke ein Glied mit kritischem Index (* - i'ii) erzeugt wird, 

+— I 



P. Stankiewioz' Buchdruokerei, Berlin Sjfcy IIHfJUJwa^erstr. 14. 
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